Microscopic model for the higher-order nonlinearity in optical filaments 
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Using an exactly soluble one-dimensional atomic model we explore the idea that the recently 
observed high-order nonlinearity in optical filaments is due to virtual transitions involving the con- 
tinuum states. We show that the model's behavior is qualitatively comparable with the experimen- 
tally observed cross-over from self-focusing to de-focusing at high intensities, and only occurs at 
intensities which result in significant ionization. Based on these observations, we conjecture that 
this continuum electron nonlinear refraction exhibits strong memory effects, and most importantly, 
the change of its sign is effectively masked by the de-focusing due to free electrons. 



Since the first experimental observation of long- 
distance propagation and optical filamentation in high- 
power femtosecond light pulses more than a decade 
ago [H, it has been accepted that the main nonlinear 
effects controlling the phenomenon are the optical Kerr 
effect and de-focusing due to the free electrons generated 
by high-intensity ionization. However, Loriot et al. re- 
cently presented an experimental measurement of higher- 
order, intensity-dependent Kerr nonlinearity in optical 
filaments @, y], which exhibits a cross-over from self- 
focusing to de-focusing at high intensities. This develop- 
ment was quickly followed by a theoretical work Q which 
concluded that the long-standing theory of optical fila- 
mentation in gases needs to be changed in radical ways. 
In particular, it was proposed that the occurrence of free 
electrons is not necessary for formation of femtosecond 
filaments. At time of this writing, the experiment still 
awaits an independent corroboration. Moreover, a mi- 
croscopic explanation of the higher-order nature of the 
nonlinear refraction beyond the usual Kerr effect has not 
been offered so far. Using full quantum simulations of an 
atom subjected to a short pulse Nurhuda and co-workers 
have attempted to deduce the form of the higher-order 
nonlinear refraction (see e.g. Q). However, these efforts 
are hampered by the difficulty of separating the various 
contributions to the total nonlinear optical response. 

In this paper we advance the idea that higher-order 
nonlinear refraction can arise from virtual transitions 
from the ground atomic state to continuum states and 
back to the ground state. We term this continuum elec- 
tron nonlinear refraction in contrast to the more usual 
bound electron nonlinear refraction that involves virtual 
transitions from the ground state to bound states and 
back to the ground state jfjj]. The bound electron non- 
linear refraction is usually calculated within third-order 
perturbation theory and leads to the familiar Kerr non- 
linearity for transparent dielectrics. In contrast, as we 
shall see, evaluation of the continuum electron nonlinear 
refraction involves a non-perturbative calculation that 
leads to a nonlinear optical response that can cross over 
from self-focusing to self-defocusing at high intensities. 

To elucidate the physics of continuum electron non- 



linear refraction we employ an exactly soluble one- 
dimensional atomic model where the electron-ion inter- 
action is modeled using an attractive delta-function po- 
tential, sometimes termed 'delta-Hydrogen'. This model 
has the virtues that in the absence of any external fields 
it has only one bound state, the ground state, plus the 
continuum states, and the model remains soluble even 
in the the presence of an external field. This allows 
us to isolate the continuum electron nonlinear refraction 
since there is no bound electron nonlinear refraction in 
this model. We remark that the one-dimensional atomic 
model has been explored extensivel y by the mathematics 
0, H| and physics communities (9l4l3|. Previous works 
involving light-matter interactions include tunneling ion- 
ization [3l in strong laser fields [IU [l6| and also the 
time dependence of the survival probability of the decay- 
ing ground state |17j . 

The occurrence of higher-order nonlinear refraction has 
clear wide ranging implications for the dynamics of self- 
focusing collapse and the formation of filamentation in 
general, and so a goal of our study is to pose and answer, 
at least partially, the following questions raised by the 
Loriot 's experiment: 

A) What is the microscopic mechanism behind the in- 
tensity dependent nonlinearity? 

B) Which occurs earlier, ionization or the higher-order 
Kerr focusing-defocusing cross-over? 

C) To what extent is the higher-order nonlinearity in- 
stantaneous? 

D) Perhaps most importantly, why no previous experi- 
ments produced a clear evidence of the higher-order 
nonlinearity? 

In this work, we address the above problems in the 
low-frequency limit, i.e. when the ratio of the photon 
energy to the ionization potential is sufficiently small so 
that we can study the metastable ground-state in a quasi- 
static approximation. While it is obvious that the simpli- 
fied one-dimensional atomic model employed here is too 
simple to produce any quantitatively verifiable outputs, 
we contend that it provides insight into the underlying 
physics of higher-order nonlinear refraction. This is sup- 
ported by the fact that higher-order nonlinear refraction 
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appears to be universal, and in particular common to 
atoms and molecules, and as such it should only depend 
of the most basic properties of the given system. 

For the sake of reader's convenience, we next recall 
the definition of the model, and its exact solution in the 
form of the Hamiltonian's resolvent. Then we identify 
the relevant resolvent pole in the non-physical sheet of 
the spectral parameter, and from it we calculate the field 
dependent polarizibility of the ground state. 

In a symbolic form the Hamiltonian can be written 
as is usual in the physics literature, namely in terms of 
a delta-function potential acting on a one-dimensional 
particle in a homogeneous external field: 



H = 



-h 2 d 2 



-h 2 A 



S(x) — eFx 



(1) 



2m dx 2 2m 
More precisely, the domain of the Hamiltonian consists of 
locally absolutely continuous functions, ACi oc (K), with 
their derivative in AGi oc (IR\0), which satisfy: 



dip(x) 



dx 



dip(x) 



dx 



-Aip(0) 



(2) 



a:=0- 



x=0+ 

Here, the quantity A, which we assume to be positive, 
measures the strength of the delta-function potential. 
The above condition "encodes" the delta-function in the 
domain of the Hamiltonian whose action is then defined 
through the differential operator 

-h 2 d 2 i/>(x) 



Hil){x) = 



eFxij){x) 



(3) 



2m dx 2 

which is the same as in the system with no contact poten- 
tial. It is also required that the result of §5$ is quadrati- 
cally integrable. To simplify notation, we utilize units in 
which m = 1/2, e = 1, and h = 1. 

The spectral properties of H are fully captured in its 
resolvent (H — A) -1 or, equivalently, in the Green's func- 
tion G{x,x',X) = (x\(H - A)"V) 



G[x, y, A) = G Q (x, y, A) 



G o (x,0,X)G Q (0,y,X) 



(4) 



1/A-G (0,0,A) 

where Go(x, x' , A) is the Green's function for A — 0. The 
above is an exact result and a consequence of the Krein's 
theorem Q. We refer the reader to [l6| for an intuitive 
derivation. 

With A = 0, Go (a;, a;', A) satisfies the equation 

- d xx G (x, y, A) - (Fx + X)G (x, y, A) = S(x - y) , (5) 

and can be expressed through a pair of solutions to the 
homogeneous equation: 

r (a- „ \\ - -^L(x,X)tp R (y,X)/W(t/j L ,ip R ) , x < y 
Uo{x,y,A)- _ MxA)My)X y w ^ LM 5 x>y 

where the Wronskian W(^ L ,ip R ) = i/j L d x ^ R - ip R d x i/) L , 
and the left and right solutions iI>l,r must tend to zero 
at their respective infinities 

lim %j) L (x,X) = Q, lim %l> R {x,X)=0, Im{A} ^ . 
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Re[A] 



FIG. 1: (Color online) Contour plot of zeros of the real (black) 
and the imaginary (red) parts of Resolvent resonances ap- 
pear where the contours intersect. The row of solutions close 
to the real axis correspond to the zeros of Ai. Poles along 
e -2 !7 r/3 are c i oge ^ zeros f Ci + . The resonance that origi- 
nates in the bound state is marked by the circle in the right 
panel. On the left, the field is weak, and the resonance lo- 
cation is barely distinguishable from the ground-state energy 
E g = -1/4. Left: F = 0.025, right: F = 0.125. A = 1. 



In the upper half-plane of the spectral parameter, 
Im{A} > 0, we can write ipL,R in terms of Airy func- 
tions as follows: 

il> L (x, A) = Ai(-0 x ) = Bi(-0 + iAi(-0 

with (assuming F > 0) 

£ = F 1/3 (x + X/F) and W = -F 1/3 /ir 

The asymptotic expansions appropriate for the respective 
sectors of the complex plane 18( show that at infinity 
these solutions behave as 



ip L - exp (-2^ 3/2 /3) and ip R ~ exp (+2^73) 

and therefore vanish exponentially for Im{A} > and x 
approaching negative and positive infinity, respectively. 
Thus, tpL.R is indeed the appropriate pair of solutions 
to calculate the Green's function. Equivalently, they 
can be viewed as combinations of the Hankel functions 
H^ 2 3 (2/3£ 3/2 ). This may be advantageous for "travers- 
ing" the entirety of the two spectral-parameter sheets 
since all needed analytic continuations can be conve- 
niently obtained with a single pair of functions. How- 
ever, we will restrict our attention to a sector in which 
the above Airy representations are sufficient and easy to 
use. 

Having obtained explicit expressions needed for the full 
resolvent (j4| we now turn our attention to its poles. The 
denominator of the second part gives the equation to find 
the resonances, namely 1 = ^4Go(0, 0, A) which explicitly 
reads 



3/2, 



nA „ 
1 = — tttAi 



^1/3 I p2/3 



-X 



;Vi{ fA) 



Bi 



I F 2 / 3 



(6) 
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This equation has infinitely many solutions; for a small 
F, zeros of Ai and Ci + = Bi + iAi give rise to two fami- 
lies of resonances. However, we are interested in another 
solution which converges, with F — > 0, to the zero-field 
ground-state energy E g = —A 2 /4 on the real axis while 
approaching from the lower half-plane. Figure [TJ illus- 
trates all three types of solutions in weak and strong 
fields. The resonance corresponding to the metastablc 
ground state is marked by a circle. 

The real and imaginary parts of this resonance pole 
location Ar = E r — iE\ are related to the ionization rate 
and nonlinear polarizibility, respectively. Specifically, the 
ionization rate is proportional to its imaginary part 



1 

w = — 

T 



2Ei 

h 



A useful scaled quantity is the ionization rate in units of 
time given by the oscillation period of the driving field. 
If we specify the wavelength equivalently as the number 
n of photon energies needed to overcome the system's 
ionization potential \E g \, then the corresponding scaled 
ionization rate reads 



U'r, 



2n- 



Fk 
\E„ 



This quantity tells us if we should expect an appreciable 
fraction of atoms to be ionized within a single cycle of 
the optical field. For argon with its ionization potential 
of 15.7 eV and for the wavelength of 800 nm, n = 11. We 
will therefore use w\\ to place our results in the context 
of femtosecond filamentation. 

To relate E r to the nonlinear polarization, we express 
the ground-state energy shift E r — E g as the energy of an 
induced dipole in the electric field. In turn, the induced 
dipole is expressed through the ficld-dcpcndcnt polariz- 
ibility a(E): 



E, 



E g =pE = a(E)E 2 



The nonlinear component of the polarizibility, namely 
«nl(-E) = a(E) — a(0) is the quantity we aim to compare 
to the results of the higher-order Kerr measurement. As 
with the ionization rate, it is convenient to utilize a scaled 
form of the nonlinear polarizibility, namely 

Xnl = a- NL (E)/a(0) 

This quantity, termed scaled nonlinear susceptibility, 
measures the induced susceptibility of refraction relative 
to the linear-regime susceptibility of the unperturbed sys- 
tem. 

We next discuss the properties of the exact solution to 
the resonance equation ([6|). We assume that the sys- 
tem has an ionization potential of 15.7 eV, i.e. that 
of argon, the incident pulse wavelength is 800 nm, and 
we express the field strength F equivalently through the 




intensity [TW/cm 



FIG. 2: (Color online) Nonlinear susceptibility (black) and 
ionization rate (red) obtained from the exact solution for the 
resolvent pole. 



light intensity /. Figure [5] depicts the result obtained 
from the exact solution. The nonlinear susceptibility 
exhibits nearly linear increase for low intensity, which 
corresponds to the usual Kerr nonlinearity when the in- 
duced index of refraction is proportional to the light in- 
tensity. At higher intensities, xnl(-0 saturates and sub- 
sequently decreases into negative values. The zero cross- 
ing occurs at I s» 45TW/cm 2 . For very high intensi- 
ties (not shown) the scaled susceptibility approaches — 1 
which corresponds to complete cancellation of the atom's 
linear susceptibility. However, the extreme high-intensity 
regime is practically irrelevant because the survival prob- 
ability of the ground state vanishes rapidly. The second 
curve in Fig. [2] shows the scaled ionization rate for the 
800 nm driving wavelength. We can infer that in the 
vicinity of the susceptibility zero crossing, already a ten- 
cycle pulse is sufficient to reach ionization probability of 
about eighty percent. We can also see that when the 
nonlinear susceptibility starts to saturate, ionization al- 
ready sets in. In other words, the nontrivial behavior of 
the nonlinear susceptibility is "synchronized" with sig- 
nificant ionization. 

For a complementary view of the interplay between 
ionization and polarizibility, it is instructive to exam- 
ine asymptotic solutions of ([6]) for weak fields. Previous 
works [H| [Til showed that the ionization rate is non- 
perturbative, and has the functional form of the tunnel- 
ing ionization. This means that there exist no Taylor 
expansion in terms of the field strength. Because the 
polarizibility and ionization are connected through the 
resolvent pole location, it is only natural to expect that 
also the nonlinear index will be non-perturbative and lack 
proper Taylor expansion; this is indeed the case as illus- 
trated next. In the lowest order (for A = 1) the pole 
equation and its resonance solution Ar become [l6[ 



1 = 



e sf 



Ar . = - 



1 



2V\ 4VA 4 
The imaginary part gives rise to the non-perturbative 
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tunneling ionization rate 14J. Here one can see that 
there are actually two small parameters, namely F and 
e _ 6F. We develop the asymptotic solution of © in 
powers of these. For fixed expansion orders, we obtain 
a converged approximation in a finite number of itera- 
tions. The procedure shows that the real part of Ar will 
also contain non-perturbative terms, starting with e~w. 
This means that neither ionization nor the polarizibil- 
ity can be developed as a Taylor expansion in the field. 
A second important point is that even very high-order 
asymptotic solution can only reproduce the polarizibility 
curve roughly up to 10TW/cm 2 in Fig.[2j This is because 
the "small" parameter e - ^ becomes large very quickly 
with increasing field. This sheds light on the fact that 
the measured nonlinear coefficients seem to represent a 
"divergent" function. The representation in powers of 
intensity, ^ri2fc/ fc , is actually ill suited to represent this 
nonlinear behavior, and in our opinion contributes to the 
rather large error bars. 

Finally, let us summarize our finding from the exact 
solution of the model system, and relate them to the 
questions posed in the introduction. 

A) Most importantly, we can see that even a model 
with a single bound state exhibits an intensity depen- 
dent nonlinearity very much similar to that observed by 
Loriot et al., and that it occurs in the comparable in- 
tensity range. The absence of other bound states in the 
model indicates that the effect is due to virtual transi- 
tions into continuum. We speculate that, similarly to 
the high-harmonic generation, the few necessary "ingre- 
dients" are limited to the existence of ground state and 
of a continuum spectrum with the band-gap setting the 
scale for the intensity. 

B) We showed that the nontrivial behavior in the po- 
larizibility is intimately connected with the ionization, 
which sets in as soon as the susceptibility starts to sat- 
urate. Moreover, in this respect it is important to real- 
ize that in a real system the excited bound states will 
give rise to the additional n^I Kerr nonlinearity. This 
will then shift the real part of the total susceptibility to- 
ward higher intensities, while leaving the imaginary part 
intact. Consequently, the plasma generation will start 
before the nonlinearity changes its sign to de-focusing. 
From here we conjecture that in naturally occurring fil- 
aments, the higher-order nonlinearity will be effectively 
masked by free electrons. We speculate that this may 
very well be the answer to our question D), namely why 
previous experiments did not noticed these effects. 

C) The fact that it is the continuum states that give 
rise to the nontrivial nonlinear behavior strongly suggests 
that the high-order nonlinearity will exhibit finite mem- 
ory whose influence is likely to become stronger with in- 
creasing intensity. This is because unlike discrete states, 
the continuum of states can form "sufficiently many" su- 



perpositions which can "encode and remember" the his- 
tory of the system. Of course, to verify this conjecture, 
time-resolved investigations will be necessary. If corrob- 
orated, a departure from a truly instantaneous (or very 
fast) nonlinearity will have an important consequence in 
decreasing the conversion efficiency in generation of high 
frequencies and supercontinuum. 

To conclude, while keeping in mind the simplicity of 
the studied model, we believe that it gives us very im- 
portant clues that will be eventually useful for under- 
standing microscopic mechanisms controlling filamenta- 
tion processes on few-femtosecond time scales. 
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